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We calculate local potential energy minima (inherent structures) for a simple model of orthoter- 
phenyl (OTP), from computer simulations over a wide temperature range. We show that the results 
are very sensitive to the system size. We locate, from the temperature and size dependence of in- 
herent structure energies, a well defined cross-over temperature that corresponds to the dynamical 
critical temperature predicted by mode coupling theory and in mean field spin glass models. Com- 
paring the known phase space structure of mean field p-spin models and the present results for the 
OTP model, we find evidence supporting the proposed similarity between glass-forming liquids and 
a class of spin glass models. 
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The idea that essential features of glass formation from 
supercooled hquids [Q,^ may be contained in the behav- 
ior of spin glass models has been actively pursued in 
recent years, with a focus on mean field p-spin models 
. In particular, the dynamical ideal mode coupling 
theory (MCT) for supercooled liquids |^ and the disor- 
dered mean field p-spin models, with p > 2 display many 
common features observed in structural glasses. Both 
approaches predict, on cooling the system, the existence 
of a dynamical critical temperature (Tmct in MCT lan- 
guage and Td in p-spin language) at which, from a dy- 
namical point of view, an ergodic to non-ergodic tran- 
sition takes place. Below this temperature the system 
is dynamically frozen in a metastable state and confined 
in a finite region of phase space (basin, in the follow- 
ing). Relaxation to the true equilibrium state below this 
temperature takes place via activated processes, which 
are not taken into account dynamically in MCT and are 
absent in the mean field p-spin case, where all basins 
are separated by infinite barriers. Both theories make 
predictions for the time dependence of correlation func- 
tions (correlators). The evolution of correlation functions 
in the p-spin case coincides P,p[ with that obtained for 
the schematic MCT model [BL In real structural glasses 
no sharp ergodic to non-ergodic transition is observed, 
due to the finite height of the barriers separating differ- 
ent basins, so that this dynamical critical temperature 
plays the role of crossover temperature between regions 
of different dynamical character. Specifically, it has been 
argued that below this crossover temperature, activated 
crossing of potential energy barriers, first suggested by 
Goldstein [|l0|, becomes the dominant mechanism of dif- 
fusion MCT has been the subject of several tests 
above Tmct and it has been shown to describe the space 



and time dependence of correlators satisfactorily ||T2|,|l3| . 

On the other hand, mean field p-spin models have the 
merit of offering a clear picture of the structure of phase 
space, in particular the aspects associated with the tran- 
sition at Tjj. The analytic solutions of the TAP equa- 
tions 1^,0 show that the region of phase space visited 
below T]j is composed by an exponentially large number 
of different basins, separated by infinitely high barriers. 
Different basins are unambiguously specified by the value 
of the energy at the minimum, Emin, which assume val- 
ues only between a minimum value Eo and a maximum 
value Eth- Stationary solutions with Emin higher than 
Eth are unstable even at zero temperature. Below To, in 
thermodynamic equilibrium, basins with a well defined 
value of Emin{T) are populated at each temperature. 

The TAP solutions aXT — K coincide with the local 
phase space minima of the potential energy and are thus 
analog of the so-called inherent structures proposed by 
Stillinger [|l5| to describe the phase space structure of liq- 
uids. In Stillinger's approach, phase space is partitioned 
in basins such that steepest descent trajectories starting 
from any point in a basin converge to the same energy 
minimum, which defines the inherent structure config- 
uration. Previous numerical studies of inherent struc- 
tures in supercooled liquids have been confined to rather 
small system sizes and have rarely aimed at identifying 
the presence of an underlying dynamical critical temper- 
ature. A recent study of a Lennard- Jones binary mixture 
at constant density for a 256 atom sample did not 
reveal any particular signature of the presence of a dy- 
namical temperature in the T dependence of the inherent 
structure energies. 

This Letter explores the possibility that the phase 
space structure of the liquid state is similar to that of p- 
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spin models. The work is based on a detailed study of the 
energy relaxation in phase space toward the minimum, 
starting from equilibrated independent configurations at 
a series of temperatures. The use of the conjugate gradi- 
ent technique [ pT[ guarantees an efficient monotonically 
decreasing path toward the minimum and of course the 
absence of basin change during minimization. The ba- 
sic idea behind this work is that the equivalence between 
inherent structures of the liquid and TAP solutions fa- 
cilitates the comparison of the phase space structure of 
structural glasses and of mean field p-spin models. 

We perform a series of long molecular dynamics (MD) 
simulations for systems of different sizes ranging from 
343 molecules (the smallest size compatible with the in- 
termolecular potential range) up to 9261 molecules [^ . 
The selected intermolecular potential models the OTP 
molecule as a rigid structure composed of three sites. All 
site-site interactions are described by the same Lennard- 
Jones potential, whose parameters have been optimized 
to reproduce the properties of the liquid. For this poten- 
tial Tmct — 280 K, a figure calculated from the T depen- 
dence of the non-ergodicity parameter [|9j. Below Tmct, 
the MD trajectories have been calculated for times up 
to 20 ns, to guarantee thermodynamic equilibrium and 
satisfactory sampling of phase space. Both static and 
dynamic quantities, evaluated along the MD trajectory, 
do not show any finite size effect. Several different equi- 
librium configurations, separated by more than the re- 
laxation time have been used as starting points for the 
evaluation of the corresponding inherent structures by 
performing a standard conjugate gradient minimization. 
The energy minimization has been performed up to a rel- 
ative precision of 10^^^ to ensure a proper determination 
of the minimum. 

Fig. shows the potential energy during the mini- 
mization procedure at different temperatures. While for 
T < 305 K there is one single relaxation process to the 
minimum, with a characteristic "time" which has no sys- 
tem size dependence, at T = 346 K the relaxation to the 
minimum occurs in a two-step process. The first (size 
independent) relaxation process is followed by another 
slower relaxation which drives the system to a deeper 
and size dependent minimum. These data indicate that 
the system first relaxes quickly to a quasi-stable size in- 
dependent point characterized by an energy Es{T) and 
then explores the possibility of finding a lower energy 
minimum configuration with value Emm{T). Strikingly 
enough, the success of this search is strongly size depen- 
dent. Indeed, for the smallest system EgiT) coincides 
with Emin{T). It is particularly important also to note 
that the change of behavior from size dependent to size 
independent relaxation to Emin occurs very close to the 
previously estimated Tmct- Inp-spin models, the energy 
Emin of the disordered solution (i.e. the stable solution 
for T > Td) is equal to that of the first (highest) sta- 
ble TAP solution. By extrapolating Emin{T) vs 1/A^ 



for \/N going to zero (see inset of Fig. |] we find that 
limjv^oo Erain{T = MQK) is close to Emin{T = Tmct) 
which thus provides strong evidence for the identification 
of Td with Tmct- 

In p-spin models, below To the energies of the relevant 
TAP minima lie below the threshold value Eth, and the 
curvature in phase space of these TAP solutions is pos- 
itive in all directions. Stationary solutions with energies 
greater than Etu display unstable (or non-positive curva- 
ture) directions |20 ^. If this feature is preserved in the 
phase structure of supercooled liquids, one would predict 
that, in the thermodynamic limit, for T > Tmct there 
will be always a finite number of unstable directions for 
higher energy stationary points and during a local en- 
ergy minimization, the system will always relax to Eth- 
But if the fraction of unstable directions /„ is sufficiently 
small, as it is close to Tmct, in finite size systems there 
would be realizations of higher energy states which are 
stable. For finite size systems, only when 6Nfu{T) is 
greater than one, i.e. only when statistically there will 
be at least one unstable solution, the system can aban- 
don the saddle state and relax toward a lower energy 
minimum. To test the validity of this prediction. Fig. 
H shows the energy during the minimization procedure 
at a T close to Tmct for all studied realizations of the 
large size system. For comparison, we show also the re- 
laxation curves averaged over all realizations for the 343 
molecules cases. Interestingly enough, a clear bimodal 
effect is observed. A fraction of the realizations relaxes 
to Emin, while the remaining fraction is trapped at the 
energy Es- The smaller system on average relaxes to 
Es- As suggested by Fig. |^, such system size depen- 
dence is seen only above the previously estimated Tmct 
for this system. The size effects displayed in both fig- 
ures, as well as the dichotomous behavior observed in 
Fig.g, suggests that Tmct can be identified as the tem- 
perature at which the fraction of unstable directions in 
phase space /«(T = Tmct), evaluated at Eg, goes to zero 
from above. In Ref. [pT|, by analysing the instantaneous 
normal mode spectrum ||22| in a deep supercooled liquid 
above Tmct, it was shown that the fraction of truly un- 
stable directions in phase space tends to vanish at Tmct, 
strongly suggesting that the dynamics above and below 
Tmct is dominated by different physical transport mech- 
anisms: Above Tmct the dynamics is controlled by the 
search of unstable directions while below it is controlled 
by activated processes. We believe that our present ob- 
servations are relevant in establishing the nature of such 
a dynamical crossover. 

Fig. ^ shows the temperature dependence of Es{T) 
and Emin{T) for different system sizes. Es{T) does not 
show any significant dependence on the system size, sug- 
gesting that in searching for the minimum the system 
quickly finds a local phase space environment that does 
not depend on the system size. For T higher than a 
size dependent cross-over temperature, a secondary re- 
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laxation process appears - suggesting that statistically 
at least one direction in phase space at Es{T) is unsta- 
ble, such that the system relaxes to the deeper energy 
minimum Emin{T). The cross over temperature is size 
dependent and in the limit 1/N going to zero extrap- 
olates to T ~ 2SbK (see inset of Fig. ||), supporting 
the identification, in the thermodynamic limit, of this 
cross over temperature with Tmct- Above the cross-over 
temperature there is a large T range where E,nin{T) is 
approximatively constant, suggesting that there is a well 
defined minimum characteristic energy of the liquid state 
above Tmct- Below the crossover temperature the differ- 
ence between Es{T) and Ejnin{T) disappears. No finite 
size effects are present. The phase space visited in this 
T range is partitioned in many different basins. Interest- 
ingly enough, the variance in Emin{T) is very small below 
the crossover temperature, supporting the view that in 
the thermodynamic limit, at each temperature, basins 
with a well defined value of Erain{T) are populated. It 
is also worth noting that Es{T) is not only system size 
independent, but is an increasing function of T and coin- 
cides with Ejriin up to the crossover temperature for each 
system size. The size effects shown in Fig. ^ suggest the 
study of the phase space of the liquid not only in terms of 
minima, but also in terms of stationary points of higher 
order, as phase space features relevant to the system's dy- 
namics. Finite dimensional cuts of the thermodynamic 
phase space (finite N) transform saddle points into lo- 
cal minima, altering the description of phase space as a 
more numerous collection of different basins. In the ther- 
modynamic limit, only the region of phase space visited 
for T < Tmct can be described as a collection of differ- 
ent basins, in agreement with the description of the free 
energy of supercooled liquids proposed by Stillinger p^ ] 
and observed in the TAP solution of the p-spin models. 
Following the analogy with disordered p-spin models, we 
note that Emin{T) correspond to the equilibrium TAP 
solutions. Interestingly, the finite size effect effects we 
report allow visualization of the equivalent of the TAP 
solutions for E > Eth, which have been shown to be un- 
stable stationary states |2C[] . 

In summary, in this Letter we have presented results on 
the phase space structure of supercooled liquids (struc- 
tural glasses) and discussed the analogies with disordered 
spin glasses. In particular (i) we have identified Tmct 
in terms of a static property of the system, showing that 
in the thermodynamic limit, Tmct corresponds to the 
temperature at which the T dependence of the inher- 
ent structure energy has a clear break. This suggests a 
possible experiment based on ultrafast quenches to de- 
termine the Tmct for a liquid, (ii) we have shown evi- 
dence that above Tmct the inherent structure energy is 
T- independent , suggesting the possibility of a unique 
basin characterizing the liquid, (iii) Above Tmct, size ef- 
fects become important due to the transformation of low 
dimensional saddle points to minima in finite systems. 



The different structure of the phase space visited above 
and below Tmct demonstrated by the present work offers 
precise hints about the difference of dynamical behavior 
above and below Tmct- While below Tmct the slow- 
ing down of dynamics on cooling can be connected to 
the time required to hop across energy barriers between 
basins, above Tmct exploration of phase space can take 
place via unstable directions connecting different quasi- 
stable environments (defined by their value Es{T)). This 
suggests further investigation of such quasi-stable solu- 
tions for their relevance to the dynamics above Tmct 
. (iv) Finally, based on the formal equivalence between 
TAP solutions at T = QK and inherent structures, we 
have shown that the phase space structure of the present 
model liquid is consistent with the phase structure of the 
mean field p-spin models and have identified T]\jct with 
Td- 
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Indeed, it appears plausible that the onset of stretched 
exponential and non-Arrhenius behavior above Tmct 
studied in is related to Eg (T) , which is under current 
investigation. 



FIG. 1. Relaxation of the potential energy per molecule during the minimization procedure for three different system sizes, 
starting from equilibrated configurations at three different temperatures T — 266, T = 305 and T — 346 K. Each curve is the 
average over several independent minimizations. Dotted line. A'' = 343; Full line A'^ — 1000; Long Dashed N = 9261. The inset 
shows the inherent structure energy as a function of the inverse of the number of molecules N for T = 305K (dashed line) and 
for T = 3467^ (fuU line). 



FIG. 2. Potential energy per molecule as a function of the number of steps in the conjgate gradient procedure when 
N = 9261. Different curves refer to different equilibrium configurations at T = 318A". The inset shows an enlarged view in the 
energy region around Es (full lines) together with the averaged relaxation for N = 343 (long dashed). 



FIG. 3. Temperature dependence of Es (empty symbols) and Emin (filled symbols) as a function of temperature for three 
different system sizes. Square N — 9261; Circle N = 1000, Triangle N — 343. Lines ar drawn to guide the eye. It is important 
to note that (i) below 300 K, Es coincides with Emin', (ii) at T = 3Q5K finite size effects hinder the second relaxation process 
also when N = 9261. 
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FIG. 1. F. Sciortino et al 
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FIG. 2. F. Sciortino et al 
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FIG. 3. F. Sciortino et al 
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